Dynamic instability transitions in ID driven diffusive flow with nonlocal hopping 
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One-dimensional directed driven stochastic flow with competing nonlocal and local hopping events 
has an instability threshold from a populated phase into an empty-road (ER) phase. We implement 
this in the context of the asymmetric exclusion process. The nonlocal skids promote strong clustering 
in the stationary populated phase. Such clusters drive the dynamic phase transition and determine 
its scaling properties. We numerically establish that the instability transition into the ER phase is 
second-order in the regime where the entry point reservoir controls the current and first order in the 
regime where the bulk is in control. The first order transition originates from a turn-about of the 
cluster drift velocity. At the critical line, the current remains analytic, the road density vanishes 
linearly, and fluctuations scale as uncorrelated noise. A self- consistent cluster dynamics analysis 
explains why these scaling properties remain that simple. 

PACS numbers: 64.60.Cn, 05.70.Ln, 05.40.-a, 02.50.Ga 



I. INTRODUCTION 

Driven stochastic diffusive flow of particles in nar- 
row channels is widely used as a prototype process to 
study dynamic scaling. The central interests in such pro- 
cesses reside with the structure of non-equilibrium sta- 
tionary (NES) states, the dynamical pathways to those 
NES states, and the scaling properties of dynamic phase 
transitions. Driven stochastic processes often undergo 
dynamic phase transitions inside their NES states as a 
function of control parameters. This is true even in one 
dimension where equilibrium phase transitions are for- 
bidden. One-dimensional (ID) processes serve thus as 
ideal platforms to unearth novel general principles af- 
fecting the dynamic scaling properties and the structure 
of NES state distributions. The maximum entropy prin- 
ciple familiar from thermal equilibrium, leading to the 
Gibbs distribution, does not apply to driven systems. We 
need to formulate new principles, if possible, to predict 
the structures of NES states ^ . Such principles can 
emerge from the study of a wide array of specific pro- 
cesses, in particular using numerical simulations and/or 
exact solutions. 

In this paper, we address how clustering induced by 
nonlocal stick-slip events affects ID driven stochastic dif- 
fusive flow. Our model is a generalization of the con- 
ventional asymmetric exclusion process (ASEP) I, i. 
We allow nonlocal forward hopping events across empty 
stretches of road, which introduces an instability towards 
an 'empty-road' (ER) phase. 

Clustering and queuing has been studied recently in 
different generalizations of ASEP. In the conventional 
ASEP, particles hop stochastically in a preferred direc- 
tion along a loop, without being able to pass each other 
and with the occupation of each site limited to one. That 
ASEP stationary state is surprisingly trivial, completely 
disordered. However, it is very sensitive to defects and 
boundary conditions. The sensitivity to reservoirs in 



open boundary type set-ups was discovered first. In that 
set-up the stationary state is known exactly by the so- 
called matrix product ansatz 0, 0] . The two reservoirs at 
the edges of the channel compete with the bulk for control 
over the bulk density p and the average flow rate. This 
gives rise to phase transitions between bulk and reservoir 
controlled phases. The latter can be viewed as elemen- 
tary forms of queuing (traffic jams). The density profiles 
near reservoirs have exponential or power-law tails. The 
scaling properties remain rather simple and are quali- 
tatively understood from the fact that the ASEP bulk 
stationary state is fully uncorrelated, while fluctuations 
spread as / ~ t^/^ with z = 3/2, and travel with group 
velocity Vg = 1 — 2p (see e.g., Ref. @). 

Introducing bulk impurities to ASEP, e.g., a single 
slow-bond in the bulk, represents the next level of com- 
plexity in queuing. Such a single static defect creates 
more intricate types of jamming than edges, because, 
unlike reservoirs, information and correlations can travel 
across the obstacle. Recently we demonstrated that the 
scaling properties of the slow-bond queuing transition are 
non-trivial. We established the existence of an interme- 
diate phase with a pure power-law shaped queue, and 
showed that the transition into the macroscopic jammed 
phase does not take place until a finite strength of the 
obstruction |6|, [3]. 

A next step in the queuing and clustering saga is to 
allow such impurities to be mobile and participate in 
the dynamics. Macroscopic clustering represents jam- 
ming. For example, in the two species process by Arndt 
et al [1] particles of opposite charge hop in opposite di- 
rections and the passing of species creates jamming. The 
stationary states of such processes show strong cluster- 
ing. The possibility of phase separation into macroscopic 
clustered states is discussed at length in the literature 
both numerically [9| and analytically^^ as well as the 
link to the zero-range process (ZRP) 

In this paper, we study queuing caused by a different 
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mechanism. In our generalization of ASEP, clustering 
is induced by nonlocal hopping events. A particle can 
jump with probability 1 — p to the nearest-neighbor site 
(if empty) or with probability p all the way forward to the 
site immediately behind the particle in front of it. We use 
open boundary conditions with reservoirs at both edges. 
Starting-off as uniform at p = 0, the stationary state 
becomes increasingly lumpy with well-defined clusters as 
p increases. Then, at a critical value Pc, this clustered 
liquid becomes unstable and the road empties out. The 
location and the nature of this phase transition are our 
main concerns. The transition turns out to be discontin- 
uous or critical. Clustering and coarsening prove to be 
essential in explaining both. 

Clustering is quite common in driven stochastic flow. 
It appeared for example in the two species type process 
by Arndt et at. Q. Those clusters were successfully de- 
scribed by the ZRP [Til, which have zero drift. In our 
process, the nonlocality induces clustering with nonzero 
drift. More importantly, that aspect induces first- and 
second-order phase transitions, and the scaling proper- 
ties of those are governed by the clustering. 

This process and its phase transitions also relate to 
various viscous and/or dissipative stick-slip type phe- 
nomena, where the particle number is the only locally 
conserved variable. Potential applications include traffic 
and/or granular flow [12], the motion of stick and flow in 
sandpiles [13], phase separation in steady sedimentation 
of colloidal crystals 0, ll5|, electronic and/or molecular 
transport in nanoscale systems [l^, phase ordering in 
rough films [13], the motion of molecular motors driven 
by ATP [18^, the motion of a depinned flux lattice in a 
current-carrying superconductor [19] and so on. Nonlocal 
hopping can be used to mimic the dynamic features asso- 
ciated with, e.g., the competition between maximum and 
minimum speed/drift in traffic/granular flow, inertia of 
falling grains inside avalanche processes in sandpiles, and 
the role of gravity in dynamic and static sedimentation. 

This paper is organized as follows. In section [TTl we 
present numerical results for the phase diagram, clus- 
tering, and the nature of the phase transitions. In sec- 
tion \H1\ we setup a self-consistent free cluster analysis, 
in which the stationary state is described in terms of a 
collection of freely drifting clusters, i.e., self-organized 
mesoscopic collective objects that absorb and emit indi- 
vidual particles constantly. The first-order segment of 
the phase transition line is explained as the result of a 
turn-about in the drift velocity of those free clusters. The 
critical segment of the transition line is caused by star- 
vation and the fluctuations in the density near the entry 
reservoir become crucial. To describe this, we extend 
in section IIVI the self-consistent cluster analysis to the 
mother cluster, which is the cluster attached to the en- 
try point reservoir. This analysis explains the numerical 
details of the critical transition presented in section 
and discussed in section [Vll The paper concludes, in sec- 
tion IVIIi with a brief summary of the results and some 
open questions. 
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FIG. 1: (Color online) a—p phase diagram from our numerics. 



Our self-consistent cluster approximations resemble re- 
cent "zero-range process" (ZRP) treatments of clustering 
phenomena, e.g., in the ASEP with two or more species 
of particles (driven in opposite directions) [10]. Our ap- 
proach is both more basic and more general. In those 
processes the clusters are stationary on average, while 
ours have a non-zero drift velocity. In the ZRP approach 
the internal structure of the clusters is ignored, while in 
our approach the internal density profile is taken into ac- 
count. Both aspects are essential for the transition into 
the ER phase. ZRP can be solved exactly, while our 
discussion remains more approximative. 



II. PHASE DIAGRAM 

Consider a chain of length 1 < x < L with sites that 
can only be empty or occupied, Ux = 0, 1. The chain is 
open, in contact with the entry side reservoir at x = 
and the exit side reservoir at x = L + 1. The updated 
rule for our process is as follows: 

(1) Select a site, < x < L, at random. 

(2) If X and it is occupied, the particle slides to 
its nearest neighbor site x-\-l (with probability 1 — p) or 
detaches and jumps all the way to the empty site directly 
behind the nearest particle in front of it (with probability 

p). 

(3) In case the chain is completely empty in front of it, 
the particle jumps all the way forward into the exit side 
reservoir. 

(4) If the entry side reservoir, x = 0, is selected and 
site X = 1 is empty, a particle jumps onto site x = 1 
with probability a. Nonlocal hopping events from the 
reservoir are not allowed in our set-up. 

Figure [1] shows the phase diagram obtained from our 
Monte Carlo simulations with road lengths up to L = 
16 000. At p = the process reduces to the conventional 
ASEP. There, the stationary state is well-known to be 
in either: the input reservoir-controlled phase (C) for 
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FIG. 2: (Color online) Spatiotemporal density profiles for L — 256 (1 < x < running horizontally from left to right) up to 
t = 512 (0 < t < running vertically from top to bottom) illustrate clustering in both (a) (a = 0.25 : pc — 0.19) the C and (b) 
{alpha = 1.0 : Pc — 0.30) the MC phase. Initial conditions include one compact free cluster in the middle of the system to 
show clearly that the drift velocity of the free clusters turns around at p :^ 0.3. 



a < ^, where the bulk density is set and limited by the 
influx fugacity a; or the bulk-controlled phase (MC) for 
q; > ^, where the influx rate is not the limiting factor 
anymore and the bulk density is set by maximizing the 
bulk current. 

The C and MC phases persist into p > 0, but be- 
come increasingly clustered, see Fig. [2l Their histori- 
cal names become quite unfortunate. 'MC phase' stands 
in the conventional ASEP context, for 'maximal current 
phase'. Control over the current, by either the bulk or 
the reservoir, is closely linked at p = with the so-called 
'maximal current principle'. The latter has been pro- 
posed as a general principle to predict dynamic station- 
ary states. Its success in the conventional ASEP is closely 
tied to the stationary state being uniform and clusterless, 
and well- approximated by mean-field (MF) theory. The 
maximal current principle arose naturally as a reformu- 
lation of the MF iteration process. Its generalization to 
strongly clustered states remains unclear, and this issue 
plays an important role in the discussion below. 

'Bulk-in- Control' and 'Reservoir-Controlled' would be 
better names for the two phases in Fig. [21 but we avoid 
confusion by preserving their historical abbreviations 
'MC ('the-Middle-in-Control' phase) and 'C ('being- 
Controlled'-phase). In the C phase, the stationary-state 
current varies both with p and a, while in the MC phase 



it only varies with p regardless of boundary conditions, 
i.e., a-independent; see Fig.[3l This feature distinguishes 
the two phases in our numerical analysis. 

For strong nonlocal hopping (large p), the stationary 
state transforms into an 'empty road' (ER) phase. The 
transition line starts at small a as a critical line, located 
just below p = a. It levels-off for increasing a, and 
changes at the critical endpoint, at a = 0.64(1), into 
a first-order transition line. The latter is a strictly hor- 
izontal line, with pc = 0.300(2), see Fig. [H The density 
of particles pr jumps to zero at the first-order line, see 
Fig. m but the average current remains continuous, see 
Figs. [3] and El 

The road is not truly empty in the ER phase. The 
average bulk density vanishes, but typically a finite clus- 
ter of particles 'hangs' still from the entry side reservoir 
near x = 1, and some isolated clusters and individual 
particles are traveling through the bulk. Figs. [6] and [71 il- 
lustrate this. They show the evolution of the probability 
distribution for finding Nr = prL particles on the road 
across the first-order transition, at a = 0.75, and across 
the critical line, at a = 0.2. Note the gradual evolution 
of the distribution across at a = 0.2 versus the abrupt 
change at a = 0.75. 

Consider the average current through a bulk bond be- 
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FIG. 3: (Color online) The steady- state current J versus p 
for L = 8000 at fixed a = 0.1, 0.15, • • • ,1 from bottom to top 
(with the same a-labelling as in Fig.H)). 

tween x and x -\- 1 

= (n^^x+i) -^p{vxVx^i) -pPo{x + 1), (1) 

and through the edge bonds 

(Ji) = a{vi) , (J^+i) = {ul) ^p{vl) -pPo{L), (2) 

where Ux is the occupation operator at site x, and Vx = 
1 — Tlx is the vacancy operator. Equation ([1]) states that 
the current through the bulk bond is zero if site x + 1 is 
occupied; equal to one if site x + 1 is empty while x is 
occupied; and equal to p if both sites are empty provided 
at least one particle exists on the road somewhere to the 
left of X + 1. The latter requires the introduction of the 
vacancy string operator 

X 

Po{x) = {\[vy), (3) 

which counts the probability for the entire road to the left 
of X + 1 being empty. Po{L) acts as an order parameter, 
which vanishes in the C and MC phases, but remains 
non-zero across the entire chain in the ER phase. Po(^) 
is always non-zero near the entrance at small x. 

The qualitative structure of the phase diagram is now 
easily understood. In the stationary state, the average 
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FIG. 4: (Color online) The steady-state bulk road particle 
density pr versus p at fixed a = 0.1,0.15, ••• ,1.0 for L = 
8000. 
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FIG. 5: (Color online) J — p versus p at fixed a = 
0.2, 0.25, • • • , 1.0 for L = 8000. The intersection points with 
the J — p = line represent the C-ER and the MC-ER phase 
transition points. 

current is constant throughout the system; and Eq. ([2]) 
implies the identity 

a{v,) = l-{l-p){vL)-pPo{L) (4) 

This excludes the C and MC phases from extending into 
the p > a region of the phase diagram. The identity 
cannot be satisfied for p > a with Pq{L) = 0; the right 
hand side varies between p and 1, because < < 1, 
and that is inconsistent with < 'Ui < 1 for p > a. 

Equation ([2j) also implies that the current along the 
entire transition line is equal to J = Pc{ol)] because at 
Pc-, {vl) is already equal to one, while Po(^) is still equal 
to zero. Our numerical simulations confirm this. 

The levelling-off of the transition line into a horizontal 
(constant p) segment is linked to this as well. As a start, 
we consider MF theory. In the MC phase, where the 
bulk controls the current, Eq. ^ is then approximated 
as J = Vh — {l—p)v1 where Vh is the vacancy density deep 
inside the bulk. The maximum current representation of 
MF theory, SJ/Sv^ = 0, yields, Vb = ^ and J = with 
r = 1 — p. This together with the condition J{pc) = Pc 
puts the MC-ER transition line at p = ^. 

This is only an upper limit to its true location. The 
clustering in the MC phase, shown in Fig. [21 is not rep- 
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FIG. 6: (Color onUne) The probability distribution P{Nr) 
for finding Nr particles on the road, at a — 0.75 for various 
p = 0.295, 0.296, 0.303 (from right to left) at L = 8000. 
The MC-ER transition occurs at pc ^ 0.299. 
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FIG. 7: (Color online) The probability distribution P{Nr) 
for finding Nr particles on the road, at a = 0.2 for various 
p = 0.154, 0.156, 0.164 (from right to left) at L = 8000. 
The C-ER transition occurs at pc ~ 0.161. 

resented in MF theory. It underestimates the number of 
nearest neighbor pairs. These pairs shift the transition 
line downwards, see Eq. ([T]). More advanced versions of 
MF theory incorporate local correlations but still fail to 
account for clustering. Instead, we develop in the follow- 
ing section a self-consistent cluster approach. It predicts 
the correct location of the MC-ER transition line, and 
more importantly, it provides insight in the mechanism 
that makes the MC-ER transition discontinuous. 

We note that a closed chain (with periodic boundary 
conditions) does not include the ER phase due to con- 
servation of the total number of particles and holes in 
such a set-up. A similar model was studied in the con- 
text of the mass-conserving coalescence process, the so- 
called mass-chipping model [U, [22[ with and without a 
bias. Implications of our results in this context will be 
discussed elsewhere [23|. 

III. FREE CLUSTER ANALYSIS 

Adopt the point of view where the C and MC type sta- 
tionary state configurations are fully clustered; with one 
distinct 'mother cluster' attached to the entry reservoir 
at small x and many disconnected freely drifting clusters 
in the bulk. These clusters are meant to be only meso- 
scopic in size, because true phase separation in terms of 
macroscopic high density (clusters) and low density sec- 
tions is not realized in this process, see Fig. [2l 

Assume the bulk clusters are compact non-fractal 
mesoscopic objects with average bulk density pc and lo- 
cal average densities, pp and pR, at the front xf and rear 
xr end, respectively. The front moves forward by creep 
(local hops, Xf ^ xf and backward by detachment 
of the front particle by a nonlocal hop to the next cluster, 
with velocity 

" = ~ ^) ~ p/^^ 

The rear moves by local forward hops (internal nonlocal 
hops are suppressed by the high particle density inside 



the clusters) and by attachment of a new particle due to 
a nonlocal hop from the free cluster behind it. 

UR = — = [1- PR) -p (6) 

The drift velocity Ud of the cluster is associated with its 
center of mass 

Ud = ^{uf ^ur) = l-p - ^pR - ^p/pF- (7) 

A cluster maintains its integrity only if its internal bulk 
density pc is stationary. The number of particles in the 
cluster, A/'c, is determined by nonlocal hops at its two 
ends, and is thus on average constant in time, dNc/dt 
—p -\- p = 0. A stationary internal density pc requires 
therefore that the average cluster length (^=xf — XRhe 
invariant. Equations ([5]) and (|6]) then relates the front 
and rear densities as 

PR- PF = P- (8) 

Inserting reasonable values for the densities yield imme- 
diately that Ud changes sign as p increases, and does so 
before the MF estimate of Pc. For example, pF = 1/2 im- 
plies that the drift velocity changes sign at Pc = 1/3; and 
also that pR = 2/3, large enough to validate the cluster 
concept self-consistently. 

This turn-about of the drift velocity of the clusters ex- 
plains the first-order nature of the MC-ER transition. 
Any initial configuration decays (after forming clusters) 
to an empty road for Ud < 0^ because all bulk clusters 
travel backward and fall back into the entry point reser- 
voir. This change in drift velocity is clearly visible in the 
density profile time evolutions of Fig. [2l 

Emboldened by this success, we dare to push the self- 
consistent cluster analysis further. We need to relate 
the above equations to the free cluster current Jc and 
bulk density pc to obtain a closed set of equations. To 
determine the free cluster properties throughout the C 
and MC phases, and also estimate the value of Pc- 

Visualize each cluster as a separate conventional 
ASEP, i.e., the absence of significant numbers of internal 
nonlocal hopping events due to the high local internal 
density. In conventional ASEP the density probability 
distribution deep inside the bulk is spatially uncorrelated. 
It is safe to presume the same is true inside the clusters. 
This yields 

Jc = (1 - rvc)vc, (9) 

with r = 1 — p. The cluster acts like an individual ASEP 
with open boundary conditions in a somewhat unusual 
set-up: It has a fluctuating 'lattice size' ^ and a fixed 
injected current instead of a fixed lattice size and a 
fluctuating injection current (controlled by a fugacity like 
a). The total current inside a cluster is equal to the 
injection current plus the contribution of the center of 
mass drift as 

Jc=P^PcUd- (10) 
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Jc and pc should not be confused with the global current 
and density. They only coincide with those at the tran- 
sition point Ud = Eq. (p!Q|) reproduces correctly the 
exact result that J = Pc at the transition. 

Equations ([7j)- (p!Q|) provide still only 4 conditions be- 
tween 5 variables: the cluster current Jc, the drift veloc- 
ity Ud, and the rear, bulk, and front densities {pR^ pc^ 
and pf)- Typically (at a given value for p), a set of pos- 
sible solutions can be found. Let's parameterize them by 
their values for Jc. We use the maximal current princi- 
ple to close the equations and select the correct value of 
Jc from the set. The justification of this is again based 
on conventional ASEP, where the exact matrix method 
solution as well as the MP approximation have this prop- 
erty [3]. A secondary feature of the conventional ASEP 
is also preserved, that the stationary state density profile 
is always uniformly decreasing, pn > Pc ^ Pf- 

This procedure leads to two distinct type free cluster 
solutions as function ofp: a 'front-limited' (F) state (with 
pc = Pf) Sit p > Ps] and a 'maximum-current' (MC) state 
(with SJ/Svc = 0) dit p < Ps. The crossover takes place 
at Ps = 1(3 -Vb) =0.191. 

The MC free cluster solution at p < resembles very 
closely and smoothly connects to the p = global MC 
state. It has internal vacancy density '^c = ^5 cluster 



current Jc = j-, drift velocity Ud 



|, and edge 



densities pp = 2p and pR = \. However, we cannot 
take the MC free cluster solution seriously at small p, 
because the internal cluster density becomes too small 
for the cluster concept to remain meaningful. Instead, 
we interpret ps as marking the onset of crossover towards 
the non-clustered MC stationary state at p = 0. 

The F solution with pc = Pf, maximizes the internal 
current inside the cluster for p > pg. The cluster density 
is equal to 



Pc 



1 



2(1 -p) 



V4p— — p 



(11) 



The rear density, pR = p/pc^ is larger than pc = pF 
in the entire interval Ps < p < Pci reflecting a uniformly 
decreasing density profile (required for internal stability). 

These P solutions tell us that inside the MC phase 
(bulk-in-control phase) of Pig. [H near Pc, the configura- 
tions are clustered and that each of these free clusters 
is limited and controlled by its front. Inserting p^ = pp 
from Eq. (dH) into Eq. ([5]) yields Ud = Oditpc = 1/3. This 
estimate for the first-order transition line is remarkably 
close to our numerical results, pc — 0.30. 



IV. MOTHER CLUSTER ANALYSIS 

We now turn our attention to the the critical line be- 
tween the C and ER phases. This transition is dom- 
inated by the fiuct nations near x = 1. Therefore, we 
extend the self-consistent cluster analysis to the mother 
cluster, i.e., the cluster attached to the entry point reser- 
voir. This analysis explains the scaling properties of the 



second-order transition, as determined numerically and 
presented in the following sections. 

The mother cluster is presumed to be a mesoscopic 
object, connected to the entry side reservoir. It governs 
the bulk current in the C phase, while being controlled 
itself by the reservoir. The mother cluster acts like sl p = 
ASEP with a fiuctuating front, i.e., with a fiuctuating 
length. The injection current from the reservoir in its 
rear is equal to Jc = avi^ with Vx = I — px before. 
The bulk density of the mother cluster can be presumed 
to be uncorrelated again, Eq.([9]). The evolution of its 
front dit xf obeys Eq. ([5]). An amount p of the current 
through the cluster fiows away at the front of the cluster 
by nonlocal hops, and the remainder extends the cluster 
forward: 



Jc=P^ PcUf = p + Pc [(1 - p) - p/pf] . 



(12) 



This set of equations is still incomplete. It ignores the 
probability Po{x) that all sites to the left of a specific site 
X are empty, see Eq. ([1]). In the free cluster analysis this 
played no role {Po{x) = 0) because another particle can 
always be found to the left of the free cluster as long as 
the mother cluster exists and remains mesoscopic. The 
mother cluster current is very sensitive to the density 
profile near the reservoir, and thus to Po{x). 

The mother cluster density profile can be presumed 
(self-consistently) to be sufliciently uniform that we can 
apply mean field theory internally. 



Vx+i = Jc/[l-rvx -pPo{x) 
Po{x) = v,Po{x-l). 



(13) 
(14) 



Compare this with Eq. ([T]). MP theory for conventional 
ASEP gives erroneous power-law and exponential density 
profile tails. Similarly, the following analysis can be at 
best only qualitatively correct. 

Solutions are found by solving the above equations it- 
eratively and by finding the value of Jc that satisfies the 
boundary conditions at both the rear and front. It is 
important to realize that the solution is unique for 'fi- 
nite' cluster lengths ^c- The mother clusters are finite 
in length, but we analyze them as if they are mesoscopic 
and effectively infinitely long. Por infinite cluster lengths, 
the equations typically allow a range of Jc with possible 
solutions. The maximum current solution is the correct 
one, because it coincides with the ^c ^ limit of the 
unique finite length solution. 

At ^c = CO, the construction of a solution involves the 
matching of two iteration processes: the forward itera- 
tion, X ^ X + 1, starting from the rear of the cluster and 
the backward iteration, x x — 1, starting from the front 
of the cluster. It is useful to sketch briefiy some of the de- 
tails of this. The iteration process has two bulk densities 
as possible fixed points, determined by the roots of the 
quadratic equation Jc = Vc{l — rVc) (using that Pq ^ 
in the bulk of the mother cluster). The low density fixed 

point Vc"^ is unstable and the high density fixed point Vc^^ 
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FIG. 8: (Color online) Forward iteration processes. 

is stable in the forward iteration process. They reverse 
roles in the backward iteration process. 

For a given value of Jc, the forward iteration process 
starts at site x = 1 with density vi = Jc/o^, and typ- 
ically iterates towards vi^-' (see Fig. [8j). The backward 
iteration process starts at site x = xf with density vf 
prescribed by Eq. ([12]), and typically converges towards 
Vc"^ . Those solutions do not match, unless the current is 
raised to the value where the two fixed points coincide. 
This matching is equivalent to (and the origin of) the 
maximal current principle. This type of solution, where 
the two fixed points merge, applies to the bulk controlled 
MC state. It exists when the starting densities of the for- 
ward and backward iteration processes are located inside 
the attraction domain of the vi^"^ and Vc"^ fixed points, 
respectively. 

This fails in the front and rear controlled, F- and R- 
type states. In the F states, the starting point of the 
backward iteration process, the value vf^ crosses over to 
the unstable side vf < of the fixed point, while 
we raise the current. In that case Jc must be chosen 
such that those two densities match, vf = vi^\ The 
density profile at the front is thus fiat. The rear profile is 
typically able to iterate forward towards this fixed point 
Vc^^ , and thus match the backward iteration. (This solu- 
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FIG. 9: Mother cluster states in the a-p phase diagram. 

tion still maximizes the current, in the boundary limited 
sense.) 

In the rear-controlled (R) states, the roles of the for- 
ward and backward iteration processes are reversed, i.e., 
the forward iteration process crosses over to the wrong 
side of the Vc"^ fixed point. This requires more care than 
the F state, because the initial value vi = J/a needs to 
be iterated forward into the region x Ip where Po(^) 
becomes negligible, before it becomes clear where this 
profile 'lands' compared to the low density fixed point 
Vc^ . The latter is the reason why we need to bother with 
the iteration process for the mother cluster, while for the 
free clusters it suffices to write down the boundary and 
bulk conditions, and use the maximal current principle 
formulation of the matching process. 

Figure [9] shows the mother cluster states as a func- 
tion of a and p. The MC (bulk-controlled) and F (front- 
controlled) states are basically identical to the free cluster 
states discussed in the earlier section. This confirms the 
stability of the MC clustered phase at large a. In these 
states the front or bulk of the mother cluster is in control, 
setting Jc, and at the rear the reservoir is able to supply 
the requested current without letting the mother clus- 
ter detach. The F state can be divided into two states, 
F-up and F-down, distinguishing between whether the 
forward iteration process approaches from below or 
above (see Figs. [8] and [9]). In the free cluster analysis 
only F-down is present. Due to the presence of Pq{x)^ 
the mother cluster is not required to be uniformly de- 
creasing any more. 

These F and MC solutions describe mesoscopic mother 
clusters. The growth velocity (drift velocity of the front) 
iii? is positive. Strictly speaking the mother cluster could 
keep growing until it fills the entire chain, but that would 
lead to the uniform MF-type non-clustered state along 
the entire chain. In reality, density fiuctuations make 
the mother cluster break apart regularly, emitting free 
clusters. 

The transition into the R states at small a takes place 
when the front (F-state) or bulk (MC-state) imposed 
value of Jc can not be satisfied at small x anymore. From 



there on, the rear takes over as the hmiting factor, and 
the front end of the profile is required to follow instead. 
This is the origin of the starvation process responsible for 
the second-order transition line. 

The mother cluster front can only follow the rear at 
small p, denoted in Fig. [9] as the R-up and R-down states. 
The current Jc must be tuned such that the forward it- 
eration process lands exactly onto the low density fixed 
point Vc^. In the R-up and R-down states, the back- 
ward iteration process, starts at a 'Ui? inside the attrac- 
tion basin of the low density fixed point, and therefore 
the backward iteration process is able to converge. 

The R-up state has an increasing density profile in the 
front, where Pq = 0, and is therefore dynamically un- 
stable. Moreover, the bulk densities in both the R-up 
and R-down states are very small, just as in the MC 
state, and self-inconsistent with the cluster concept. The 
mother cluster density profile starts high at x = 1, but 
decays to the small unstable fixed point value beyond 
X = Ip. Such clusters must be unstable to internal clus- 
tering fluctuations and break-up regularly into a small Ip 
sized object. 

Everywhere near the second-order transition line, in 
the area marked as R in Fig. [9l the backward iteration 
process fails because vp starts at a < vi^^ outside 
the attraction basin of Vc\ The forward and backward 
iterations can not be matched anymore. A true R-type 
mother cluster solution seizes to exist. The front of the 
mother cluster can not be in a locally stable stationary 
type forward creeping state anymore. 

This inability to construct a stable R cluster front near 
the second-order transition line in Fig. [H leaves no other 
options than to toss out the mother cluster front con- 
ditions, Eq. (p!2|) . and redefine the mother cluster as a 
quasi- stationary R-type object (typically short, except 
very close to the second order transition line) defined by 
the half-space forward iteration process only, with Jc the 
precise value of the cluster current that let vi = a/Jc 

iterate into Vc\ The mother cluster size is set by the 
length scale Ip where Po{x) and the density becomes 
small. Ip is now the onset of a turbulent low density re- 
gion where the formation of new free clusters takes place, 
continuously without any "local stationary state order" 
or stability. 

Next, we face the problem that we lost our (maximum 
current) justification for why Jc must settle on that spe- 
cial value where the forward iteration process converges 
into the unstable fixed point. Nevertheless, it is arguably 
the dynamically stable choice in the following sense. Con- 
sider the half-space solutions with smaller values of Jc. 
Their forward iteration process converges into the sta- 
ble high density fixed point. Those clusters would have 
a pronounced dip in their density profile near Ip (when 

the iteration process sails near the Vc'' fixed point) be- 
fore increasing towards the large density value of the vi^^ 
fixed point. Such structures are unstable and tend to 
shrink. They resemble mother clusters in the process of 
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FIG. 10: (Color online) The road density pr versus p at a = 
0.2 for various L. The inset illustrates the finite-size behavior 
of pr at the C-ER transition point in a log- log plot. 

pinching-off a free cluster. The breaking- up event shrinks 
the mother cluster towards Ip. Moreover, the break-up 
also increases the current of the combined two objects. 
Half-space solutions with values of Jc larger than the spe- 
cial choice, make the iteration process crash almost im- 
mediately through zero density, and can be interpreted as 
representing very short mother clusters. Those are also 
unstable and tend to grow. The empty space in front of 
such a short cluster can carry away still only a current of 
order p, while Jc has increased. Moreover, it is less suc- 
cessful in breaking- up because Po{x) is not small at its 
front, and therefore the free cluster immediately in front 
of it is subject to a fluctuating injection current, larger 
on average than p, and thus with a significantly smaller 
or even negative drift velocity. 

The C-MC phase boundary is represented in the 
mother cluster analysis as the threshold, see Fig. [9l be- 
tween the R states and the MC and F states. These 
reproduce the true location (from our numerical simula- 
tions) not too well. It has a kink instead of being smooth. 
The reason for this is likely that the mother cluster is 
very short in the R-states and badly approximated by 
the mesoscopic description. 

The C-ER critical line is represented in the mother 
cluster analysis as the line where the R-state solution 
iterates to zero density. This estimate virtually coincides 
with the true location determined from our numerical 
simulations. This remarkable accuracy is likely linked 
to the fact that the mother cluster length, of order Ip, 
diverges at the transition, and the mesoscopic approach 
becomes valid. 

The crossover from the horizontal first-order segment 
of the transition line to the curved second-order part 
represents a critical endpoint. It is located at exactly 
Pc = 1/3 and a = 2/3 in this analysis, and is also close 
to the true location. Our numerical simulations locate it 
at Pc — 0.300(2) and a :^ 0.64(1). The scaling proper- 
ties of the critical endpoint appear numerically not spe- 
cial, and consistent with the notion, emerging from the 
cluster analysis, that the origins of the first- and second- 
order transitions (starvation versus drift-velocity turn- 
about) are distinct and unrelated. Their scaling behav- 
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FIG. 11: (Color online) Details of J — p versus p as a function 
of p at a = 0.2 for various L show that finite-size corrections 
are very small and J{pc) — Pc — 0.16125(5). Compare with 
Fig. El 



iors are thus simply superimposed onto each other at the 
critical endpoint. 

Within our numerical accuracy this critical endpoint 
coincides with the endpoint of the C-MC boundary. 
There is no a priori reason to expect this coincidence, 
but in the mother-cluster approach it arises naturally. 



V. NUMERICAL RESULTS AT CRITICALITY 

In our numerical simulations, we average over 10^ ~ 
10^ independent runs up to t = 10^ Monte Carlo steps 
for non-stationary state features. For the stationary state 
properties we simulate system sizes L = 500x2"^ with n = 
1, 2 • • • ,5 and discard all data before t :^ L^, or more, 
to allow the system to reach the stationary state. We 
collect data up to t = 10^ ~ 10^ Monte Carlo steps and 
average over 10 ^ 200 independent runs. We monitor: 
the road density pr = (rir) {rir = Nr/L); the density 
fluctuations, X2 = L{{n'^) — (nr)^); the Binder cumulant 
1/4 = 1 — (n^)/(3(n^)^); and the current J. The results 
are shown in Figs. [TOl and [TTl 

The scaling properties of the C-ER transition are de- 
termined numerically, along several cuts though the tran- 
sition line. It suffices to present here the details along one 
representative example, i.e., only the a = 0.2 line. The 
critical point Pc along this line can be determined by stan- 
dard methods, from the so-called crossing points or the 
maximum points in X2 as well as the crossing points in 
U4. From these we estimate Pc = 0.161(1). As mentioned 
before, the current, see Fig. [TTl must be equal J = p at 
the transition line. This is very well satisfied compared to 
the above estimate for Pc. Turning this around, requiring 
J{Pc) = Pci improves the estimate to Pc = 0.16125(5). 

The road density pr vanishes on approach of the tran- 
sition point Pc from below in the thermodynamic limit, 
as Pr ~ with A = pc — p. The density fluctuation X2 
diverges at Pc from below as X2 but becomes zero 

in the ER phase {p > Pc). The Binder cumulant U4 con- 
verges to 2/3 in the C phase (consistent with a Gaussian 
distribution) and to —1 in the ER phase (indicating an 



exponential distribution). U4 converges very fast, while 

X2™^ shows strong finite-size effects. 

The values of the critical exponents follow by exploit- 
ing conventional finite-size scaling (ESS) theory 

Pr = L-^/V(ALV-), 

X2 = L^/^^(ALV-), (15) 
U4 = /i(AL^/"), 

where u is the correlation length exponent and /, ^, and 
h are scaling functions. Elementary scaling laws yield 
hyperscaling relations such as (2/? -\- = I, 

U4 shows the least finite-size effects. We estimate the 
value of v by collapsing the U4 data, see Fig. [121 as 
l/u = 0.50(2). Collapsing the road density data yields, 
see Fig. [H P/^ = 0.48(2) and 1/u = 0.53(3). By plot- 
ting the maximum value of x against L, we can esti- 
mate the value of 7/z/. A simple power law fitting yields 
j/iy = 0.17(3), but the data is also consistent with a loga- 
rithmic fit, i.e., X2"^^ ^ log ^7 see Fig.[T4l Hyperscaling is 
satisfied only when this logarithmic divergence is correct, 
so we are inclined to put more weight on the logarithmic 
fitting. 

The value of the Binder cumulant U4{pc) is expected to 
be universal at the transition. Our numerical estimate, 
Ua{pc) = 0.08(3), is very close to zero, which may reflect 
the half-Gaussian like shape of P{Nr^ L) at criticality, see 
Fig. [71 

VI. MOTHER AND FREE CLUSTER 
SEPARATION 

The numerical data point to rather simple values of 
the critical exponents at the C-ER transition, (3=1 
and V = 2. The R-state mother cluster properties ex- 
plain why the scaling is that simple. The first and most 
crucial observation is the lack of any feature in the cur- 
rent J(p, a) across the critical line in our numerical data 
Figs. [5l and [TTl The global current appears to be analytic, 
and totally oblivious of the presence of the C-ER transi- 
tion, moving linearly through the critical value J = p, as 
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FIG. 12: (Color online) Scaling collapse of the Binder cumu- 
lants at a = 0.2 with pc = 0.16125 and 1/i/ = 0.50. 
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FIG. 13: (Color online) Scaling collapse of pr at a = 0.2 with 
Pc = 0.16125, (3/iy = 0.48, and l/u = 0.53. 

J -pc^ A{p-pc). 

The mother cluster sets the global average current. Its 
offspring, the free clusters, cannot influence this, unable 
to talk back being distinct objects with positive drift ve- 
locity. The mother cluster is an independently fluctuat- 
ing object, and the rest of the road is slaved to it. At 
the C-ER transition, the mother becomes unable to cre- 
ate free clusters at a sufficient rate to sustain a non-zero 
average bulk density pc, but is itself unaffected by what 
is going-on in the bulk. 

The global average current J{a^p) is analytic across 
the C-ER transition, because the mother cluster current 
is analytic, and the average bulk current is equal to it, 
because (on average) the mother cluster remains attached 
to the reservoir, and (on average) it does not grow in 
length (due to the breaking-off of free clusters). 

The linear vanishing of the bulk density with expo- 
nent = 1 is a direct consequence of this analyticity of 
J(a,p). They are related by the simple equation 

J{a,p) = p(l-pr)^{p^p[J^U^P)pr (16) 
- 

Pc 

( f) ( f) 

with Pc the internal free cluster bulk density and u^J ^ 
their drift velocity (see also Eq. (p!Q|) ). Consider a specific 
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FIG. 14: (Color online) The density fluctuations, X2, versus 
p at a = 0.2 for various L. The inset shows the finite-size 
behavior of the maximum values, X2™^ ^ semi-log plot. 



bond. The first term on the right hand side of Eq. ([16]) 
represents the current through this bond from isolated 
nonlocally hopping particles when this bond is not part 
of a free cluster. The second term is the current through 
the bond when it is part of a free cluster. The probability 

for the bond to be part of a free cluster is equal to pr. 

( f) ( f) 

The free cluster parameters pc and are analytic 
functions of p only and set self-consistently by the free 
cluster equations. All functions on the right hand side of 
Eq. (p!7|) are analytic, therefore the bulk density vanishes 
linearly, with P = 1. 

The C-ER transition is a bulk phenomenon induced by 
the mother cluster. It lacks therefore independent crit- 
ical bulk fluctuations. This explains why u = 2. The 
free cluster generating break-ups in the turbulent region 
near x = Ip behave like random uncorrelated events; i.e., 
the fluctuations in the number of particles injected by 
the mother cluster into the bulk behave like uncorrelated 
noise. In the C phase, those fluctuations are screened and 
absorbed by the presence and formation of free clusters. 
On approach of the C-ER transition, free clusters become 
rare and at the transition they vanish altogether. There, 
the random break-up fluctuations travel like a pattern 
across the chain. The time of flight is proportional to 
the system size. Therefore, the bulk density represents 
(biased) random noise averaged over a time proportional 
to L, and thus scales as pr ~ L~^l^ at the C-ER transi- 
tion. 

Two more variables in our numerical data express and 
confirm the same uncorrelated fluctuations. At the C- 
ER transition, the density profile (no figure included) de- 
cays to zero smoothly near the reservoir, as a power-law, 
Pr(^) ^ x~^^ with 6r ~ 1/2, consistent with random 
noise. In the C phase, the distribution of the number 
of particles on the road, P(A^^), is a Gaussian centered 
around the bulk density p^. Towards the C-ER transi- 
tion, the maximum shifts towards Nr ~ 0{1) and at the 
transition decays as a half Gaussian (see Fig. [7j) consis- 
tent with uncorrelated noise. In the ER phase, Pr{x) and 
P{Nr) decay exponentially. 



VII. CONCLUSIONS 

In this paper, we studied how the inclusion of nonlo- 
cal hopping events affects the stationary state of ASEP. 
Clusters develop, only mesoscopic ones. There is no 
macroscopic queuing transition. Instead, there is a phase 
transition towards an "empty-road phase" , which is first- 
order from the MC clustered phase (where the bulk is 
in control) and second-order from the C-type clustered 
phase (where the entry reservoir controls the station- 
ary state). The first-order MC-ER transition is induced 
by a reversal of the group velocity of the free clusters. 
The C-ER transition has rather simple scaling properties, 
reflecting the mother cluster properties near the entry 
reservoir and that the bulk (free clusters) do not interact 
back to the mother cluster. 
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One of the remaining issues is the nature of the dy- 
namic scahng. Genericahy, fluctuations in the stationary 
state broaden in time as / ~ t^l^ . The dynamic scal- 
ing in the clustered C and MC phases, is to be expected 
to remain in the KPZ universality class, with z = 3/2. 
Within each cluster length scale we expect z = 3/2 be- 
cause each cluster acts like a simple ASEP in its own 
right. At larger length scales, the process resembles a 
driven zero-range process with a preferred direction and 
non-zero group velocity, which should also have z = 3/2 
scaling. In our simulations, we did not focus on the dy- 
namic scaling properties of the C ad MC phases, but 
exploratory numerical results for z are consistent with 
z = 3/2. The dynamic scaling at the C-ER transition is 
expected to remain simple as well because that transition 
does not involve novel intrinsic bulk fluctuations. 

Clustering and its ramifications are surely the main 
message of our study. The conventional ASEP has a 
uniform stationary state with rather trivial scaling prop- 
erties; but is unstable towards clustering and queueing. 
Such clustered states, even when involving only meso- 



scopic clusters, communicate badly with each other, and 
have a hard time developing novel type fluctuations. 
Phase transitions, induced and controlled by those same 
clusters, therefore have typically rather simple scaling 
properties, like at the C-ER and MC-ER phase bound- 
aries in our process. 

It will be interesting to test our self-consistent clus- 
ter approximations to other processes with clustering in 
cases where the drift velocity of the clusters is non-zero; 
for example, a variant of the two species process by Arndt 
et al. [8] with different numbers of opposite moving par- 
ticles. 
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